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_ ' Within a microscopic statistical description of heavy ion collisions, we investi- 

D ■ gate the effect of the medium on the formation of light clusters. The dominant 

^^ • medium effects are self-energy corrections and Pauli blocking that produce the 

Mott effect for composite particles and enhanced reaction rates in the colli- 
sion integrals. Microscopic description of composites in the medium follows 
CN ' the Dyson equation approach combined with the cluster mean-field expan- 

1^^ ■ sion. The resulting effective few-body problem is solved within a properly 

cn ■ modified Alt-Grassberger-Sandhas formalism. The results are incorporated 

in a Boltzmann-Uehling-Uhlenbeck simulation for heavy ion collisions. The 
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<^ \ number and spectra of light charged particles emerging from a heavy ion col- 
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lision changes in a significant manner in effect of the medium modification of 
production and absorption processes. 
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I. INTRODUCTION 

The description of the dynamics of an interacting many-body system is particularly 
difficult when the quasiparticle approach reaches its limits. That may be the case when the 
residual interaction is strong enough to build up correlations. An example of such a system 
with correlations is nuclear matter in heavy-ion collisions. In collisions, the nuclear matter 
is first compressed and excited and then decompressed. At another extreme of macroscopic 
scale, the nuclear matter is formed during a supernova collapse and becomes the material 
of which a neutron star is made. 

Here, we address the simplest case of the formation of correlations in a nonequilibrium 
situation: hght-cluster production in low density nuclear matter. The production can be de- 
scribed microscopically following the set of generalized semiclassical Boltzmann equations. 
In the past, the coupled set was solved numerically for nucleon fjy, deuteron f^, triton 
/(, and helium-3 fh Wigner functions within the Boltzmann-Uehling-Uhlenbeck (BUU) ap- 
proach |I|J^ . The formation and disintegration of clusters generally takes place in reactions 
involving different particles. The rates for respective processes are utilized in the Boltzmann 
collision integrals. 

Rates for collisions, whether involving clusters or not, are central ingredients of all mod- 
ern microscopic approaches to heavy reactions such as the Boltzmann-Uehling-Uhlenbeck 
(BUU) approach [JI]-Q or the quantum molecular dynamics (QMD) |§|J^. For the status of 
those approaches see Ref. [0]. 

The formation and disintegration of deuterons are the simplest types of processes in- 
volving clusters discussed above. Notably, with the exception of electromagnetic production 
and breakup {np # 'yd), these processes involve the minimum of three elementary particles 
{NNN ^ Nd) and, correspondingly, require the solution of a three-body problem in the 
medium in their description. Complication brought in by the medium is evident even in the 
kinematics, as the medium brings in a preferred frame different from the center of mass of 
the three bodies. 

In general, the collision rates are affected by the surrounding medium. The common 
procedure is to ignore the medium dependence and to utilize free-space experimental cross 
sections. That procedure in many cases was very successful. To calculate the rates, including 
the self energy shift and the proper Pauli blocking, and to study the influence of the medium 
on different observables, a generalized Alt-Grassberger-Sandhas (AGS) equation |^ has been 
derived @-15]. The effective few-body problem in matter arises within the Green function 



method |TB[ when following the cluster mean-fleld expansion [|r^ or the Dyson equation 



approach |T8[. Besides the medium-dependent cross-sections, the Mott effect [p!9| , p0| plays 



an important role in nuclear matter. Both effects are studied in this paper in the context of 

a realistic simulation of a heavy ion reaction. As an example we choose the reaction studied 

in a recent experiment by the INDRA Collaboration, namely ^^^Xe+^^^Sn at 50 MeV/A ||21|| . 

The BUU approach will be explained in the following section. The needed three-body 



reaction input from the AGS-type equations will be discussed in Sec. [II. Details of the 



derivation of these equations within the Green function approach may be e.g. found in WB 



In Sec. I^ we shall present numerical results and in the last section we shall summarize the 



paper and give conclusions. 



II. NUCLEAR REACTIONS 

Within the microscopic statistical approach, the reacting system is described in terms of 
the Wigner distributions fx for hght particles, thus for nucleon /^v, deuteron f^, triton ft, 
and helion f^. Within the cluster mean- field approximation, these distributions follow the 
set of coupled semiclassical Boltzmann equations [Q , 

dtfx + {Kx, fx} = JCTUn, fd, ft, . . .} (1 ± fx) - ICTUn, fd, fu • • •} /x, X = N,d,t,... 

(1) 

In the above, Ux is a mean- field potential and {■, ■} denotes the Poisson bracket. The upper 
sign on the r.h.s. is for the Bose and the lower for the Fermi type fragments. Conversion 
between different species is accounted for in the colhsion integrals ]C{fx, fd, ft, ■ ■ ■}■ For 
example, the deuteron gain rate /C|''™(P, t) is 

/cf "(P,t) = Jd'kJd%d% \{kP\u\hk,)\l,^,, fd{h,t)f,{k2,t)f,{k,t) 

+ J d^k J d'kid'k2 \{kP\U\kik2)\%^^j,a fx{kut)fd{k2,t)fx{k,t) 

+ J d^k J d^kid^k2d^k3 \{kP\UQ\kik2k^)\l^^^^^ 

y<fN{h,t)fx{k2,t)fx{ks,t)fx{k,t) + ... (2) 

Here we use the abbreviations fx = (1 — fx) and /d = (1 + /d) and did not write overall 
energy-momentum conserving 5-functions. The ellipsis on the r.h.s. of (0) denote further 
possible contributions from four and more body collisions (e.g. tp ^ dd, hp ^ dd) or 
electromagnetic {np ^ '-/d). The loss rate /C|f^^(P, t) is given by 

ic'rip,t) = I d'k I d%d%d% \{k,k2k,\Uo\kP)\ij,__^,^^j, 

X fx{ki, t)fx{k2, t)fx{k^, t)fx{k, t) + ... (3) 

The quantity Uq appearing in Eqs. (^ and (^ is the in-medium break-up transition 
operator for the Nd — > NNN reaction and it is calculated using the AGS equation, as 
discussed in the next section. For three nucleons in isolation, Uq determines free-space 
break-up cross section cr^^ via the equation [^p2 



aliE) = ^ ~ [d%d^k2d% \{kP\Uo\k,k2ks)\' 

\vd -vx\ 3! J 

X 27i6{E' - E) (27r)=^(5(3) (^^ ^ ^^ ^ ^^^ ^ (4) 

where \vd — vx\ denotes the relative velocity of the incoming nucleon and deuteron. This 
equation along with detailed balance is used to replace the squared transition matrix ele- 
ments by the break-up cross section [Q. The cross sections can be extracted from data with 
theory aiding in interpolating and extrapolating. This type of procedure is widely used and 
has a good chance of being successful in the low density regime. However this procedure may 
not be sufficient for heavy ion collisions particularly in the lower energy regime. The ele- 
mentary cross section that enters into the Boltzmann equation depends on the medium, e.g. 
through blocking of intermediate states in scattering or through self-energy for intermediate 
states. This has been demonstrated for three-nucleon processes in Refs. P,p!T|. 



III. FINITE TEMPERATURE FEW-BODY EQUATIONS 

The effective few-body problem in matter emerges within the Green function method |TH 



when following the cluster expansion |]T^ or related Dyson-equation approach |jT8[. In the 



cluster-mean field expansion, a self-consistent system of equations can be derived describing 
an A particle cluster moving in a mean field produced by the equilibrium mixture of clusters 
with different particle number. The self-consistent determination of the composition of the 
medium is a very challenging task that is not solved till now. We adhere to an approxima- 
tion where the correlations in the medium outside the considered cluster are neglected so 
that the embedding nuclear matter is described by the equilibrium distribution of nucleon 
quasiparticles. 

The formalism to describe three-body correlations in nuclear matter has been discussed 
in detail elsewhere 0-0 . Here, we merely repeat some of the basic results for convenience. 



Let the Hamiltonian of the system be given by 

^ = E ^«i«i + ^ E ^2(12, 1'2') alaWav , (5) 

1 2W ^ 121'2' 

where ai denotes the Heisenberg annihilation operator of the particle with quantum numbers 
indexed by 1, including spin si and momentum ki. The free resolvent Go for an A particle 
cluster is given in the Matsubara-Fourier representation by 

Go{z) = {z- Ho)-' N = Ro{z) N . (6) 

In the above. Go, Ho, and A^ are all operators in the space of A particles. For simplicity, 
an index A indicating that fact has been dropped, but may get instituted when needed. 
In (H), z denotes the Matsubara frequencies zx = Ti\/{—i[3) + /x with A = 0, ±2, ±4, ... for 
bosons and A = ±1, ±3, . . . for fermions; /i is the chemical potential and (3 = l/ksT the 
inverse temperature. The effective in-medium Hamiltonian Ho for the quasi-particles in the 
mean field is given by 

n 7,2 n 

^o = E^ + s.-E^^ (7) 



where the energy shift Ei is 



with the Fermi function /i 



Ei = EV^2(12,12)/2, (8) 



The tilde in 12 indicates antisymmetrization. The factor N in Eq. @ is associated with 
particle statistics and it is given by 

N = hh...fA±hf2...fA (10) 



where / = ! — /. The upper sign in A^ is for an odd number of constituent fermions and the 
lower sign is for an even number of fermions or bosons. Note that the commutation holds: 
NRo = RqN. 

The full resolvent G in the Matsubara-Fourier representation may be written as 

G{z) = {z-Ho~ Vy^N = R{z)N , (11) 

where the potential V^ is a sum over two-body interactions within all possible pairs a in the 
cluster, i.e. 

^ = E^" = E^2°^2", (12) 



where V^ denotes the two-body potential in Eq. (y). Note that, in consequence of Eq. (|12|), 
the effective potential is nonhermitian, V^ ^ V, and also R{z)N ^ NR{z). This leads to 
the necessity of introducing the right and left eigenvectors as e.g. done in the context of 



the four-body problem in |]T5|. For a pair a = (12) of an A particle cluster, the effective 
potential of Eq. (|^) is 

{12\N^''Wi''^\l'2') = (/1/2 - /i/2)V2(12, 1'2'). (13) 

A useful operator is the channel resolvent Ga{z) within an A particle cluster, where only 
effective pair interaction in channel a is considered. This operator may be represented as 

Gc^iz) = {z-Ho- V^y'N 

= {^z-H^- N!^V^)-^N = R^{z)N. (14) 

Using inverses to the operators R defined in (P), (pJ]), and (p^, it is possible to derive 
formally the resolvent equations. For the sake of similarity with the isolated case, the A 
particle channel t matrix Tq, is defined by 

R^{z) = Ro{z) + Ro{z)T^{z)Ro{z) . (15) 

Upon introducing Ta{z) = A'"|*T2°(2;), Eq. (|15D yields to the well-known Feynman-Galitzki 



equation 16 



T^iz) = V,^ + V,^Roiz)N^T^iz). (16) 

We remark that equations of a similar form were derived by different authors in the 



past [23,24 



Note that the above equations are also valid for the two-particle subsystem embedded in a 
larger cluster (three, four, or more particles). As for the isolated equations, the effects of the 
remaining particles appear only in the Matsubara frequencies z (energies) of the considered 
particles within the cluster. No additional statistical factors related to the larger cluster 
arise. One should note that the statistical effects generally arise in the resolvent Go and not 
in the potential V2. However, it is possible to rewrite the respective equations introducing an 
effective potential, as shown in Eq. (|1^), and instead using unchanged resolvents. Following 
the more intuitive picture of the blocking in the propagation of the particles, described by 



the resolvents, we directly find in Eq. (0) the proper expression for the t matrix which 
enters the Boltzmann collision integrals (see e.g. |25[]). 



Derivation of the three-body in-medium equation is relatively straightforward and was 



given in Refs. ITTHT^. The AGS operator Upa{z) p for the three-particle system is defined 



within the equation 

R{z) = 5p^R^{z) + Rp{z)Up^{z)R^{z) , (17) 

and it satisfies the following AGS-type equation 

Up^{z) = 6f3^Ro{z)-' + J26p^N]T]{z)Ro{z)U^a{z) . (18) 

7 

This equation includes consistently the medium effects, statistical and self-energy shifts in 
the same way as the Feynman-Galitzki equation for the two-particle correlations. We use 
above the notation 6a/3 = 1 — Sa/s- The AGS-type equation yields the three-body transition 
operator U for a three-particle cluster as well the transition operator for a three-particle 
cluster embedded in a still larger cluster. In the latter case, the effect of the other particles 
in the cluster appears only in the Matsubara frequency (energy) z. The transition operator 
is defined in Eq. (|l^) in such a manner that no additional factors A^ accompany U when U 
is further employed. This in turn guarantees that the cluster equations are valid also if they 
are part of a larger cluster. The two-body subsystem t matrix entering in Eq. ( |T8|) is the 
same as the one given in Eq. (|1^). In that way it is possible to use all results of the few-body 
'algebra' recurrently, in particular those based on the cluster decomposition property. 

The in-medium bound state equation for an A particle cluster follows from the homoge- 
neous Lippmann-Schwinger type equation 

I^jb) = Ro{Eb)V\^b) = Ro{Eb)Y.N2V^\^b) (19) 

7 

where the sum is over all unique pairs in the cluster. As shown in Ref. |T3|, it is convenient 
to introduce the form factors for the three-body bound state: 

\Fp) = Y.'^p,N:iV^\,ljB,). (20) 

7=1 

With these, one derives the homogeneous in-medium AGS-type equation 

\F^) = JZ Saf3NMRo{Bs)\Fp). (21) 

(3=1 

Since the Green functions are evaluated in an independent particle basis, the one-, two-, 
and three-particle Green functions are decoupled in hierarchy. To solve the in-medium 
problem up to three-particle clusters, the one-, two-, and three-particle problems are solved 
consistently. The procedure generates the single particle self-energy shift, given by Eq. (^, 
the two-body input including statistical factors in Eq. (0), and, finally, the three-body 
scattering state. 

When solving the equations, technical reasons force us to adopt some reasonable approx- 
imations valid at smaller effective densities. Thus, the nucleon self-energy is calculated via 

6 



Eq. (I), but in the three-body equations we employ the effective-mass approximation, i.e. 
we use 

,..^ + S-(,0.^ + E-(W3). (22) 

For the sake of simphcity and the speed of the calculation, we moreover utilize angle-averaged 
Fermi-function factors N21 

N2{p, q, Pc.m.) = j—r^ I dcos6g dcosOp d(j)gd(f)p A^2(p, q, Pc.m.) , (23) 



where p and q are the standard Jacobi coordinates ||22| and the angles are relative to P, 



Rank-one Yamaguchi potentials [^ were employed in the calculations of the deuteron 
break-up cross-section in the coupled ^Si-'^Di and the ^Sq channels which should suffice 
at low energies. Through the dependence of the Fermi functions on the center of mass 
momentum, the integral equation for the three-body system exhibits a dependence on the 
net momentum in the medium frame, in sharp contrast to the free-space case. 

IV. RESULTS 

The coupled set of Boltzmann transport equations [|l],3, for the Wigner distribution 
functions of nucleons /at, deuterons fd, tritons ft, and helions fh, is solved for the collisions 
of heavy ions. For the clusters, treated in the quasiparticle approximation, it is essential to 
consider the Mott effect [0,|^,^]. Without considering this effect the description of data 



fails badly. 

The Mott effect occurs when a bound state wave function requires too many momenta 
components in the regime already occupied by momenta of other particles (Pauli blocked). 
In Ref. [jl| a geometrical model was introduced in order to account for the Mott effect in 
a numerical solution of the transport equations (|I|). Specifically, in the model the formation 
of a quasiparticle cluster is allowed only if the overlap of the respective isolated bound state 
wave function and the Fermi sphere given in Eq. (^ is lower then a specific cut-off value 
Fcut- For the deuteron, the formation is then allowed if |I|] 

Jd'qf{q + ^)\Hq)\'<F.ut, (24) 

where Pc.m. denotes the net momentum of the bound state in the medium frame and q denotes 
the relative momentum of the nucleons in the deuteron. The cut-off Fcut needs to be specified 
and it enters as a parameter the collision simulation. For the clusters in the Boltzmann 



equation set, we use the calculated deuteron and triton Mott momenta |T3| , p7| to restrict 
the choices for the parameter -Fcut- Note, however, that the procedure is not completely 
stringent since the methods of calculating the Mott momenta, following the geometrical 
picture and the t matrix approach, are different. The differently calculated Mott momenta 
are shown in in Fig. |1| for the temperature T = 10 MeV. No bound states are possible for 
the region below a respective curves in Fig. 0. As seen in Fig. |l], the parameterization given 
in Eq. (|2^) with a cut-off value of Fcut = 0.15 fits reasonably the microscopic calculation of 



the Mott momenta for a temperature of T = 10 MeV and a wide density range. Note, that 
the cut-off values for both triton/hehon and deuteron can be chosen to be approximately 
the same. We find only a slight variation of F^^t on the temperature (-Fcut — 0.17 for 
T = 5 MeV). To illustrate sensitivity to different Fc^t, we have plotted the Mott lines for 
Fcut = 0.10 and Fcut = 0.20 in Fig. |1|. A larger value of Fcut extends the space for more 
bound states. We shall study the influence of the Mott momentum on the results of reaction 
simulations. The direct implementation of the t matrix calculation for the Mott momenta 
into the BUU simulation requires the calculation of a local temperature that is very time 
consuming. Therefore we will make use of the geometrical model, with a properly adjusted 
global value of -Fcut? to account for the Mott effect in the numerical solution of Eq. (|T]). 

In Eq. (^ we give the relation between the isolated three-particle transition operators and 
the total deuteron break up cross section. The same equation is appropriate in defining the 
medium-dependent cross-section now in terms of the medium- dependent transition operator 
U given in Eq. (^). The total deuteron break-up cross section can be used in the collision 
integrals (0) and (0), as e.g. has been demonstrated in Ref. [jl[. For illustration, we show 
in Fig. 1^ the in-medium break-up cross section at a temperature of T = 10 MeV for the 



three-body system at rest in the nuclear medium, explored to a larger extend in [^,|TT],|14 
The solid line represents the isolated cross section that reproduces the experimental break- 
up data; other lines correspond to different nuclear matter densities. The influence of the 
in-medium cross sections on characteristic kinetic quantities such as the chemical relaxation 



time is significant as has been shown in fT^ . 



For a specific heavy ion reaction, we now investigate (i) to what extend the medium 
dependence of the deuteron break-up cross section affects the observables and, in particular, 
the deuteron production and (ii) how sensitive are the observables to the parameteriza- 
tion of the Mott momenta. To this end, we consider the central collision ^^^Xe+^^^Sn at 
E/A=50 MeV and compare our theoretical results with the experimental data of the INDRA 
collaboration |^ . 



To provide a first impression on how the use of in-medium rates in the Boltzmann collision 
integrals influences the outcome of the reaction simulation. Fig. ^ shows the total deuteron 
number (gain minus loss) as a function of elapsed time. The upper line is for the in-medium 
cross sections depending on the local density and the temperature. The lower line is for the 
isolated cross sections that reproduce the experimental scattering data. Clearly, the use of 
in-medium rates leads to a significant increase in the total number of deuterons. In both 
simulations, we have used Fcut = 0.15. Though the increase in the deuteron number may 
be significant, the theoretical value of the total deuteron number may be still too uncertain 
for a direct comparison to data. This is because the effect of the heavier A > 4 clusters is 
not yet included in our simulation. 

The influence of the different medium effects on the spectral distribution of proton, 
deuteron, triton, and helion clusters is shown in Fig. ^ Solid line shows the results of 
our calculation using the medium- dependent cross-sections in the collision integrals of the 
BUU simulation. Dashed line shows the spectra obtained using the isolated cross sections 
that reproduce the experimental data. In all cases we have included the Mott effect and 
for the solid and the dashed lines we use Fcut = 0.15. To demonstrate the sensitivity 
to different Mott momenta, the dotted hne shows results obtained using Fcut = 0.20 for 
medium- dependent rates. As explained in Ref. [mrT],|Il| , the deuteron break-up cross section 



is strongly enhanced near the threshold. As a consequence, we find an enhancement of about 
30% in the deuteron number in the energy range -Ec.m. < 50 MeV if we compare the dashed 
line (for isolated deuteron break-up cross-sections) to the solid line (for medium-dependent 
deuteron break-up cross sections). A larger absorption implies a larger production rate and 
maintaining the possible chemical equilibrium down to lower freeze-out densities. Inspecting 
the dotted line, we find, as expected, that the number of clusters increases with the rise 
in Fcut- Also, simultaneously, the spectra appear steeper. 

Figure |^ compares energy spectra of the light A < 3 clusters to INDRA data renormalized 
as in Ref . |^ . To compare to our calculations we renormalize our results in the same fashion 
as data, i.e. the areas below the respective curves are normalized to the same fixed value. 
As before the results obtained with the medium- dependent cross-sections are represented 
by a solid line; the dashed line shows the results of the coupled BUU calculation using 
the isolated deuteron break-up cross sections. In both calculations we use Fcut = 0.15. 
Overall, the renormalization seems to reduce the effect of using different elementary rates in 
the collision integrals. For the deuteron, the in-medium rates lead to a marginally steeper 
shape. Differences at higher cm. energies can be attributed to statistical fluctuations. 

Considering now the three-nucleon clusters, we find the helion/triton energy spectrum 
in a reasonable agreement with the experimental data. The only in-medium effect that we 
consider for the three-nucleon clusters is the Mott effect. A possible modification of the 
three-particle cluster formation and break-up is disregarded. The observed reduction, see 
Fig. ^, of the three-nucleon cluster production in the calculation using the medium deuteron 
break-up cross-section is a response to the enhanced deuteron production. 

Finally, in Fig. |^ we show the measured values (filled circles) for the mean energies of the 
light clusters emitted in the transverse direction. While the overall tendency is well repro- 
duced, the triton energy shows some discrepancy. The medium effects are less pronounced 
in this observable. The triton energy is only slightly reduced with an enhancement in the 
parameter -Fcut- No theoretical results are given for a-particles which are not yet included 
in the simulations. 



V. CONCLUSION 

Within the microscopic transport description of heavy-ion collision djTiamics, we have 
demonstrated the influence of medium effects on some typical observables. The dominant 
medium effects are self-energy and statistical corrections, i.e. for clusters the Mott effect, and 
the change of reaction rates that lead to faster time scales. The magnitude of the in-medium 
effects depends on the density and the energy deposited in the system. The basis for mi- 
croscopic calculations of the in-medium effects is the cluster mean-fleld expansion or Dyson 
equation approach. The effective few-body equations resulting from these approaches are 
numerically solved using well established and controllable few-body techniques. The cho- 
sen example of a heavy- ion reaction is ^^^Xe-|-^^^Sn at beam energy of 50 MeV/ A. Both 
effects, Mott and rate modiflcation, affect the considered observables in a comparable fash- 
ion. Presently, we have implemented the modiflcations of rates for the three-body break-up 
and formation only. However, through the coupling of the Boltzmann equations, the num- 
ber of three-nucleon clusters is also affected. We argue that a complete treatment requires 
a similar approach for the three- and four-nucleon clusters before a decisive comparison with 

9 



experiments can be made. First calculations for the Mott effect of the a-particle solving a 



proper four-body equation has been given in ||T5 
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FIG. 1. Deuteron and triton Mott momenta Piviott shown as a function of density p at fixed 
temperature of T = 10 MeV. The sohd hne represents results of the t matrix approach. The dashed, 
dotted and dash-dotted hnes represent the deuteron Mott momenta from the parameterization 
given in Eq. (M) for three different cut-off values Fcut • 
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FIG. 2. Neutron-deuteron break-up cross section for a three-body system at rest in the nuclear 
medium. Solid line represents isolated break-up cross section. Other lines are for different densities 
at the temperature T = 10 MeV. 
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FIG. 3. Integrated deuteron number as a function of elapsed collision time when utilizing 
in-medium (upper curve) and isolated (lower curve) rates from BUU simulation with F^vx = 0.15. 
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FIG. 4. Light charged particle spectra in the center of mass system for the reaction 

^^^Xe+^^^Sn at 50 MeV/A. Calculations with in-medium Nd reaction rates are represented by 

solid lines, while those with isolated Nd break-up cross sections are represented by dashed lines. 

In both cases Fcut = 0.15 is employed. Calculations with Fcut = 0.20 and in-medium reaction rates 

are represented by the dotted lines. 
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FIG. 5. Renormalized light charged light particle spectra in the center of mass system for 

the reaction ^■^^Xe+^'^^Sn at 50 MeV/A. The filled circles represent the data of the INDRA col- 



laboration |21|. The solid line shows the calculations with the in-medium Nd reaction rates, 
while the dashed line shows a calculation using the isolated Nd break-up cross section; both with 
Fcnt = 0.15. 
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FIG. 6. Mean transverse energy of light charged fragments in the angular range of 

-0.5 < cos Bern. < 0.5. 
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